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ABSTRACT 

We present new deep and accurate space (Advanced Camera for Surveys - Wide Field Planetary 
Camera 2 at the Hubble Space Telescope) and ground-based (Suprime-Cam at Subaru Telescope, 
Mega-Cam at Canada- France-Hawaii Telescope) photometric and astrometric data for the Local Group 
dwarf irregular IC10. We confirm the significant decrease of the young stellar population when moving 
from the center toward the outermost regions. We find that the tidal radius of IC10 is significantly 
larger than previous estimates of r t < 10'. By using the /, V-I Color Magnitude Diagram based on 
the Suprime-Cam data we detect sizable samples of red giant (RG) stars up to radial distances of 
18-23' from the galactic center. The ratio between observed star counts (Mega-Cam data) across the 
tip of the RG branch and star counts predicted by Galactic models indicate a star count excess at 
least at a 3a level up to 34-42' from the center. This finding supports the hypot hesis that the huge 
Hi cloud covering more than one degree across the galaxy is associated with IC10 (jHuchtmeier 1 11973: 
iCohen I fl979). We also provide new estimates of the total luminosity (Ly ^9xl0 7 L Q , My~-15.1 
mag) that agrees with similar estimates available in the literature. If we restrict to the regions where 
rotational velocity measurements are available (rssl3'), we find a mass-to-light ratio (~10 M Q /L Q ) 
that is at least one order of magnitude larger than previous estimates. The new estimate should be 
cautiously treated, since it is based on a minimal fraction of the body of the galaxy. 
Subject headings: galaxies: dwarf — galaxies: individual (IC10) — galaxies: stellar content — Local 
Group — stars: evolution 



1. INTRODUCTION 

Photometric investigations of the stellar populations 
in Local Group (LG) dwarf galaxies provide firm con- 
straints on cosmological parameters and the unique 
opportunity to investigate galaxy formation models 
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(|Mateo 1119981 : iTolstov et^l2009t IWvsell20Tol) . In this 
context dwarf irregulars (dls) play a key role, since we 
still lack firm empirical and theoretical constraints con- 
cerning their evolution a nd the possible transition into 
dwarf spheroidal galax ies jBekki 2008; Woo et al~ll2008l : 
iKormendy et al.ll2009f ). Although the number of dwarf 
galaxies known in the LG is rapidly growing in the last 
few years, current statistics indicat e that the dls are at 
least one quarter of LG galaxies (|McConnachie et al. I 
[2001 ISannalTalll2ll0"9h . 

Among the dls of the LG, IC10 is an interesting system, 
since it underwent strong star formation activity during 
the last half billion years and it is considered the only 
LG analog of starburst galaxies. Even though IC10 has 
been the s ubject of several investi gations ranging from 
the ra dio (iWilcots fc Milier"1ll998D to th e near-infrared 
(NIR, iVacca et al. 1 1200700 the UV dHunter I [20011: 
Richer et al. I 120011 ). and to the X-rav (|Wang et al. I 
20051 ). its structural parameters and in particular its ra- 
dial e xtent are poorly defined. iMassev fe; Armandroff I 
(1995) found that the" major axis o f IC10 is ~ T . A sim - 
ilar diameter (~ 6') was found by Uarrett et al. 1 (|2003D 
using th e isophotal radii from 2MASS NIR i mages. More 
recently Tikh onov &: Ga lazutdin ova I (j2009f ). using both 
ground-based and space images, suggested that the ex- 
tent of the thick disk along the minor axis is ftil 0.5'. It 
has also been suggested bv lDemers et al. I (|2004D . using 
asymptotic giant branch and red giant branch (RGB) 
stars, that IC10 should have a halo of ~ 30' diameter. 
On th e other hand, radio measurements bv iHuchtmeier I 
(|1979t hereinafter H79) indicated that IC10 has a huge 
envelope of neutral hydrogen extending over more than 
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1 square degree (62' x 80') across the sky. 

We are also facing a significant uncertainty in the 
total mass of IC10. By using Hi regions H79 found 
Mtot ~ 1.8xlO 9 M , assuming a dis tance of 1 Mpc and a 
Holmb erg diameter of ~ 10'. Also. IShostak fc Skillmanl 
(fl98fll hereinafter SS89), using high resolution maps of 
Hi regions, measured an inclination of 45° and a max- 
imum in the rotation curve of 30 km s _1 (42 km s _1 
deprojected) and the same Holmberg diameter (depro- 
jected angular diameter ~ 13'), from which they found 
Mtot ~ lxl0 9 M Q . More recently. Ivan den Bergh I (|2000D . 
following H79, but assuming a smaller distance (660 kpc, 
ISakai et al. Ill99l . found M tot ~ 6xlO 8 M . 

2. PHOTOMETRIC DATA 

The Hubble Space Telescope (HST) data were col- 
lected with the Advanced Camera for Surveys (ACS, 
pointings a,/3) and with the Wide Field Planetary Cam- 
era 2 (WFPC2, pointings 7,(5). D ata from pointings 
a, and 7 were already presented (j Sanna et al.l 120081 
2009). Pointing S includes 24 F555W-band and 24 
i ;l 814T / F-band images of 500 s each. This pointing is 
located ~ 3' NE of the galaxy center; it is outside 
the disk identified by Uarrett et al. I ((20031 . see the red 
elli pse in Fig. 1) and inside t h e thic k disk identified 
bv iTikhonov fc Galazutdinova I (|2009fl . The ground- 
based data were collected with the prime focus camera 
(Suprime-Cam, pointing e) on the Subaru telescope and 
with Mega-Cam on the CFHT (pointing £). Pointing e 
(see Fig. 1) is located across the galaxy center and in- 
cludes both shallow (3V, 3R, 31; 3x60 s per band) and 
deep (12V, 12x480 s; 13R, 13x360 s; 23/, 23x240 s) 
images. The pointing £ (see Fig. 1) is also located across 
the galaxy center and includes 3g (3x700 s) and 3i-band 
(3x400 s) images. 

Photometry on individu al image s was performed 
with DA0PH0T IV/ALLSTAR (jStetson 1 11 9871 ). The 786 
ground-bas ed images wer e simultaneously reduced with 
ALLFRAME ((Stetson I IT994T ) : the same applies to the 392 
space images. We ended up with a catalog including 
~ 1,200,000 stars with at least one measurement in 
two different bands. The ground-based data were trans- 
formed into the Johnson-K ron-Cousins sys t em us ing the 
standard stars provided bv lLandoltl (|1983l fl992l) to cal- 
ibrate local standards. The typical accuracy is 0.04 for 
I and 0.05 mag for the V band. Some external chips of 
the Mega-Cam include a limited number of local stan- 
dards and they were not included in the final calibrated 
catalog. To provide a homogeneous photometric cata- 
log the ACS and the WFPC2 were transformed into the 
V, I Johnson-K r on-Co usins system using prescriptions by 
iSirianni et al. I (|2005l ). The typical accuracy is of a few 
hundredths of magnitude in both V and /. The conclu- 
sions of this investigation are not affected by the precision 
of the absolute zero-points. 

3. RESULTS AND DISCUSSION 

The ground-based data cover a sky area of ~ 1° x 1° 
while the high angular resolution of HST images al- 
lowed us to perform accurate photometry in the inner- 
most crowded regions. We selected eight different re- 
gions, the region "C" is located across the galaxy cen- 
ter and includes data of pointings a(ACS@HST) and 7 
(WFPC2@HST), while the region "W" is located - 3' 



from the center and includes the data of the pointing 
6 (WFPC2@HST). Regions "SI", "S2", "S3", and "S4" 
cover the corners of the Suprime-Cam data (pointing e) 
and are located at ~ 18' from the galaxy center, while re- 
gions "Ml" and "M2" are two regions of the Mega-Cam 
data (pointing £) located at ~ 30' from the galaxy center 
in the SW and NE directions, respectively. Fig. 2 shows 
the /, V-I Color-Magnitude Diagrams (CMDs) of the se- 
lected regions. Data plotted in this figure show several 
interesting features. 

i) The photometry based on HST data is deep and very 
accurate, and indeed the CMDs (regions "C" and "W") 
reach limiting magnitudes of / ~ 25.5 — 26 and V ~ 
26.5 — 27 mag. The same outcome applies to the Subaru 
data, and indeed the CMDs reach limiting magnitudes 
of I ~ 25 and V ~ 26 mag. The CMDs based on CFHT 
data are shallower with limiting magnitudes I ~ 22.5 
and V ~ 24.5 mag. 

ii) Young MS stars (18 < I < 25.5, 1 < V-I < 1.5 
mag) show a strong radial gradient. Their number de- 
creases rapidly when moving from the center to the out- 
ermost galaxy regions. A handful of them are visible in 
region "W", while the blue objects of region "S" might 
be field galaxies. 

in) The different apparent colors of RGB stars when 
moving from the center to the outermost regions further 
support the occurrence of differential reddening. We es- 
timated the redde ning of region "W" using the same ap- 
proach adopted in iSanna et al.l ()2008l ). The ridgeline of 
the RGB in this field was adopted to estimate the redden- 
ing in the regions covered by Subaru data. The ridgeline 
of the contaminating blue field stars located in region 
"S4" was used to estimate the reddening in the regions 
covered by CFHT data. We found that the reddening is 
higher along the semi-major axis (E(B-V) — 0.78 ±0.10 
mag) and attains an almost constant value (E(B-V) = 
0.63 ± 0.10 mag) in the regions covered by the Subaru 
data external to the HST data. In the regions covered by 
the CFHT data external to the Subaru data the redden- 
ing attains either similar (E(B-V) = 0.63±0.10 mag) or 
smaller (E(B-V) = 0.40 ± 0.10 mag) values in the SW 
direction (Sanna et al. 2010, in preparation). The pres- 
ence of a significant reddening variation up to 20' from 
the center supports previous investiga tions based on the 
radia l distribution of Hi regions (SS89; IWilcots fc MilleTI 
[19981) . 

iv) Old RGB stars are ubiquitous and they can be 
easily identified in HST and Subaru CMDs. To prop- 
erly identify RGB stars we also plotted an a-enhanced 
isochrone (green line) of 13 Gyr at fixed metal con- 
tent (total metallicity, [M/H]=-0.66 d ex; helium content, 
Y=0 251) from the BaSTI data base (jPietrinferni et al. I 
20060- We also adopte d the same true di stance modu- 
lus (/i=24.60±0.15 mag. ISanna et al.ll2008l ) and individ- 
ual reddening estimates. The comparison between theory 
and observations further supports the evidence that stars 
located in the range 21 < / < 26 and 1 < V-I < 2.5 mag 
are bona fide old and intermediate-age RGB stars. The 
CMDs based on these data show also strong contamina- 
tion by field stars (see the region "SI"). 

v) The CMDs based on CFHT data show the same 
field contamination as the Subaru data, and probably a 

19 See also http://www.oa-teramo.inaf.it/BASTI 
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Fig. 1. — Location of ground-based and space pointings. The background is a reference image of IC10 based on randomly selected 
stars from the Digitized Sky Survey (DSS) catalog. The polygon labeled "C" shows the sky area covered by pointings a (ACS@HST) 
and 7 (WFPC2@HST), the polygon "W" the pointing S (WFPC2@HST). The large rectangular polygon shows the pointing e (Suprime- 
CamQSubaru), and the largest irregular polygon the pointing £ (Mega-Cam@CFHT). The boxes labeled "SI", "S2", "S3", " S4" and "Ml", 
"M2" display the regions adopted for the CMDs. The red ellipse displays the estimate of major and minor axis according to 1 jarrc tt et al. I 
2003). The blue circle has a diameter of 13'. 



small overdensity of stars in the region across the tip of 
the RGB (TRGB, 21 < I < 22, 2 < V-I < 3.5 mag). 

The above results indicate that the radial extent of 
IC10 has been significantly underestimated, and indeed 
according to the Subaru data the diameter is at least of 
the order of 36-46' and probably larger than one degree 
according to the CFHT data. 

To further constrain the radial extent of IC10, we 
decided to compare the observed star counts with star 
counts of foreground field stars predicted by Milky Way 
(MW) models. This approach presents several advan- 
tages when compared with the method based on the sta- 
tistical subtraction of an external control field, i) It is 
not affected by reddening differences between the galac- 
tic and the control field, ii) It is not affected by com- 
pleteness problems of the control field, thus saving tele- 
scope time. Hi) The real radial extent of these stellar 
systems it is not known in advance. Therefore, the con- 
trol fields might still be located inside their halo. The 
main drawback is that MW models need to be validated 
using de ep and accurate st ar counts covering broad sky 
regions (jRevle et al. I [2009. and references therein). 

However, the ground-based and space data sets are 
characterized by different limiting magnitudes. To pro- 
vide a robust estimate of the completeness of the former 
data sets, we compared their luminosity function (LF) 
with the LF of the WFPC2 data for pointing "W" . We 
adopted this approach since we are interested in esti- 



mating the completeness down to 1.5 magnitudes fainter 
than the TRGB, i.e., / < 23 mag. In this magnitude 
range the WFPC2 data are minimally affected by com- 
pleteness problems. We have found that the complete- 
ness is ~ 65% for the Subaru at the limiting magnitude 
of I = 23.0 mag, while it is - 65% for the CFHT at 
the limiting magnitude of I = 22.6 mag. We have cho- 
sen the above limiting magnitudes to apply conservative 
completeness corrections to both the Suprime-Cam and 
the Mega-Cam data. 

In order to compare theory and observations, we 
used the Pisa (iCastel l ani et al. I 120021: [Cignoni et al. I 
120061 ) and the Padova ([Girardi et al. 1120051 ) MW model 
We focused our attention on two outer regions cover- 
ing the same sky area, with the same mean reddening 
(E(B-V) = 0.63 ± 0.10 mag) and located in the NE 
direction, namely the "S4" (N(stars)~4000) and "M2" 
(N(stars)^ 3000) regions. The top panels of Fig. 3 show 
the /, V-I CMD for the Pisa (left) and the Padova (right) 
MW models for a sky area of 1° x 1° at the position of 
IC10 (I = 119°; b = —3.3°) and assuming a reddening of 
E(B-V) = 0.63 ±0.10 mag. 

The bottom panels of Fig. 3 show the ratio between 
the number of observed stars and the number of can- 
didate field stars predicted by the quoted MW models. 
Star counts were smoothed using a Gaussian kernel at 
fixed a. We also estimated the number of unresolved 
background galaxies, with reddened /-band magnitudes 
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Fig. 2. — I,V-I CMDs for eight regions located at different radial distances (see labeled values and Fig. 1 ). Blue young, main sequence stars 
decrease quite rapidly when moving from the center to the outermost regions. The old RGB stars are ubiquitous in HST and in Subaru data. 
The reddening also changes when moving across the different regions (see labeled values). The green line shows a a-enhanced isochroncs 
of 13 Gyr, at fixed chemical composition ([M/H]=-0.66 dex, Y=0.25) from the BaSTI database. The blue arrows display the location of 
contaminating field stars, while the red lines plotted in panel "W" and "S4" display the ridgeline used to determine the reddening. 



and V-I colors typical of the tip of RGB stars, using 
empi r ical galaxy counts (iFukugita et al. 1 1995; B enitez I 
I2000t ICapak et al. I [2004 iFerguson et al. 1 120041) and it 
was subtracted to the number of observed stars. To avoid 
spurious fluctuations caused by the limited sky area cov- 
ered by observations, theory and observations were nor- 
malized in the bright end (17.9<I<18.4 mag). The bot- 
tom left panels display the comparison with the Pisa MW 
model. There is evidence of IC10 stars across the TRGB 
region (7 = 21.66 ± 0.25 mag) in both the "S4" (left, 
1.59±0.13) and the "M2" (right, 1.38± 0.12) field. The 
star count excess is at 3cr level. The above evidence fur- 
ther supports the occurrence of IC10 stars in the region 
covered by the Mega-Cam data, since in the "S4" region 
we clearly identified IC10 RG stars. The bottom right 
panels of Fig. 3 display the comparison between obser- 



vations and the Padova MW model. There is once again 
a clear evidence of a star excess across the TRGB region 
in both the "S4" (left, 1.74±0.14) and the "M2" (right, 
1.53±0.14) field (~ 4cr level). The two Galactic models 
were constructed assuming similar input parameter^PI. 
The difference in the star count ratios between the two 
models is due to the different evolutionary inputs and to 
the normalization of the star counts in the Solar neigh- 

20 Galactic model input parameters: Kroupa initial mass func- 
tion not corrected for binaries; double exponential thin disk 
(h z (height)=250 pc, hn(scale length)=3000 pc, constant star for- 
mation rate [SFR] for t< 7 Gyr, Z(mean metallicity)=0.02); ex- 
ponential thick disk (h z =1000 pc, /ifl=3500 pc, constant SFR 
for 5<t<12 Gyr, Z=0.006); oblate halo with r 1 / 4 (h fl =2800 
pc, h 9 (semiaxis ratio)=0.8. constant SFR for ll<t<13 Gyr, 
Z=0.0002). 
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Fig. 3. — Top-Simulated CMDs for the field stars adopting a 
reddening of E(B — V) = 0.63 ± 0.10 mag in a field of view of 
l°x 1°, according to the Pisa (left) and to the Padova (right) MW 
model. Only a small fraction of the total number of stars is plotted. 
Bottom-Left-Star count ratios between the stars located in the 
"S4" and in the "M2" region with the Pisa MW model. There 
is evidence of IC10 stars at the position of the TRGB, i.e. / = 
21.66 ± 0.25 mag (see the vertical arrows). Bottom— Right— Same 
as the left, but the ratio is between observations and the Padova 
MW model. 



bourhood. 

The above findings indicate that the radial distribu- 
tion of IC10 old and intermediate-age stellar populations 
agrees quite well with the size of th e huge h ydrogen 
cloud detected by iHuchtmeier I (|1979f ) and by iCohenl 
(1979), and cover more than one square degree across 
the galaxy (r w 34 — 42'). This means that the stel- 
lar halo and the hydrogen cloud have, within the er- 
rors, similar r adial extents and resolve this pecu liar fea- 
ture of IC10 (|Tikhonov fc Galazutdinoval [20091 ) . More- 
over, this evidence further supports the hypothesis that 
the hydrogen cloud is a ssociated with the galaxy (stellar 
mass loss, pristine gas; IHuchtmeier I 119791 : [Cohen 1 119791 : 
IWilcots fc MTIIeTllSl ). 

To estimate the total luminosity we need to select can- 
didate IC10 stars. To describe the procedure, Fig. 4 
shows the /, V-I CMDs of the stars located inside a circle 
of 13' diameter across the galaxy center. The top panels 
display the photometry of space data (pointings a, /?, 7, 
S), while the bottom ones show ground data (pointing e, 
Suprime-Cam). The difference between el and e2 is in 
the mean reddening (see Fig. 4). For the stars located 
in the overlapping regions we use the HST photometry. 
The candidate IC10 stars were selected using different 
boxes in the aforementioned CMDs. The green box in- 



cludes young MS stars (V-I ~ 1.5 mag), the cyan box 
the intermediate-age stars (V-I ~ 3, 17.5 < I < 21.5 
mag), while the pink box includes old and intermediate- 
age RG stars (3 < V-I < 6, / ~ 22 mag). The position 
of the boxes in the four CMDs was shifted according to 
the local reddening (see Fig. 1). The limiting magnitude 
of the boxes is I = 23.0 according to the completeness 
experiment. 

The dashed and the dashed-dotted lines plotted in the 
bottom right panel of Fig. 4, show two young scaled So- 
lar abundance isochrones (jPietrinferni et ahl 120041 ) at 
fixed metallicity ([M/H]=-0.66 dex) and ages of t = 6 
and 200 Myr, while the solid red line shows the old a- 
enhanced isochrone with the same total metallicity and 
an age of t = 13 Gyr. These isochrones validate the 
position of the boxes we adopted to pinpoint the dif- 
ferent subpopulations of IC10. The same approach was 
adopted to select candidate IC10 stars located between 
the blue circle of Fig. 1 and the outermost regions of the 
e pointing (r < 23'). On the basis of these data and 
of t he recent IC10 dis tance based on the TRGB (830 
kpc, iSanna et al.ll2008l ) we estimated a total V-band lu- 
minosity of Ly ~ 9.13 x 10 7 Lq and a total magnitude 
of My— -15.11 mag. This estimate agrees within a fac- 
tor of two with similar es timates availa ble in the litera- 
ture (L v ~ 16 x 10 s L ? . iMateo Ill998t M v =-16.0 mag, 
I Richer et al. 1 120011 ) . The current estimate is a lower 
limit, since we are not including the stars located in the 
outermost regions covered by our photometry (pointing 
Q. However, the estimates available in the literature 
only cover the innermost galactic regions. The difference 
is mainly due to the fact that the current photometry al- 
lows us a robust identification of field stars (1 < V-I < 2, 
15 < / < 22 mag, see the blue arrow in the bottom left 
panel of Fig. 4). If they are even partially included, these 
objects introduce a systematic bias in the estimate of the 
total luminosity. The different assumptions concerning 
the adopted distance and reddening variation also help 
to explain the above difference. 

To estimate the mass-to- light (M/L) ratio of IC10 we 
restricted ourself to the galactic regions where rotational 
velocity measurements are available (see the blue circle 
with a diameter of 13' in Fig. 1). The luminosity inside 
this area is L v ~ 5.88 x lO 7 ^© (My = -14.63 mag). 
By using the quoted true distance and diameter together 
with the rotation v elocity based on Hi reg i ons ( SS89) 
we found, following IHuchtmeier fc Richter I (|2001l ) and 
iCasertano fc Shostakl ( 1980f l. a total mass of M tot ~ 
6.2xl0 8 M Q that agrees quite well with simila r esti- 
mates available in t he literature (IHuchtmeier I 119791 : 
Ivan den Bergl71l2lM IWoo et al. II2008D . Eventually, we 
found M/L ~ 10 Mq/Lq. Although this estimate is ham- 
pered by several empirical limitations it is at least one 
order of magnitude larger t han the value recently pro- 
vided by IWoo et al~l ((2008). The quoted authors use 
two independent methods to estimate the M/L ratios of 
LG dwarf galaxies: colors and inferred star formation 
history. They found that the median M/L ratio of dls 
based on the latter approach is slightly smaller than on 
the former one (0.7 vs 0.8, see their Table2). However, 
the difference needs to be investigated in more detail, 
particularly in view of the severe limitations affecting 
the estimates of the rotational velocity and of the total 
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Fig. 4. — CMDs in I, V-I bands for the stars located within 6.5' across the galaxy center. Top — space data of the pointings a, j3 and 
7 (left, E(B -V) = 0.78 ± 0.06 mag) and 5 (right, E(B - V) = 0.63 ± 0.09 mag). The colored boxes mark the CMD regions adopted to 
select candidate IC10 stars. The green box includes young MS stars, the cyan box the intermediate-age stars and the pink one the old 
and intermediate-age RGs. The position of the boxes was changed according to the local reddening. Bottom — ground based data of the 
pointing e (el, E(B - V) = 0.78 ± 0.10; e2, E(B - V) = 0.63 ± 0.10 mag). The blue arrow marks the position of candidate field stars 
(left). The dashed and the dashed-dotted lines (right) display two young, scaled Solar isochrones (t=6, 200 Myr), while the solid line an 
old a-enhanced isochrone (t=13 Gyr). Young and old isochrones were constructed at fixed total metallicity ([M/H]=-0.66 dex). 



luminosity over the entire body of the galaxy. 

We are facing empirical evidence that dls seem to 
show smaller M/L ratios when comp ared with dwarf e l- 
lipticals (see Fig. 1 and Table 2 in IW00 et a!~l 120081) . 
The new data will allow us to constrain whether this 
evidence might be affected by observational biases. 
Moreover, they can shed new lights on the prediction 



that dwarf galaxies might have ti dal radii significantly 
larger than empirical e stimates (jHayashi et al. 1 120031 : 
iKazantzidis et aLl 120041 ) . 



It is a real pleasure to thank the referee, Dr. I. Chilin- 
garian, for his constructive suggestions. 
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